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Abstract 

The polarized Gowdy model in terms of Ashtekar-Barbero variables is further reduced by including the 
Killing equations for plane-fronted parallel gravitational waves with parallel rays. The resulting constraint 
algebra, including one constraint derived from the Killing equations in addition to the standard ones of 
General Relativity, are shown to form a set of first-class constraints. Using earlier work by Banerjee and 
Date the constraints are expressed in terms of classical quantities that have an operator equivalent in Loop 
Quantum Gravity, making space-times with pp-waves accessible to loop quantization techniques. 
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I. INTRODUCTION 

This paper on the formulation of plane-fronted gravitational waves with parallel rays (pp- 
waves) in the connection formulation has two main motivations. First, space-times with pp-waves 
are interesting objects to test methods of quantizing gravity. Like midi-superspace models these 
space-times, being homogenous in two directions and inhomogeneous in the third direction, have 
a degree of difficulty lying between the complicated full theory of General Relativity (GR) and 
homogenous cosmological models amenable to quantization. Second, the recent abundance of 
investigations of conjectured Lorentz invariance violations, including an energy-dependent speed of 
light, deformations of special relativity, and experimental bounds on such effects highlight the lack 
of ab initio calculations of these effects from fundamental theories. Some of these conjectures are 
inspired by the granularity of space predicted by Loop Quantum Gravity (LQG) [l-6]. However 
the full theory seems to satisfy local Lorentz covariance [?], Q|. To our knowledge, there are no 
calculations in tractable models that clearly reconcile thegranularity of spatial geometry with 
local Lorentz symmetry. (Although Bojowald and Hossain [9] provide a derivation of gravitational 
wave dispersion in LQG in a cosmological context using perturbative methods.) Gravitational 
pp-waves appear to be simple enough for deducing a definitive answer - at least in the context of 
a model system - to the question of whether we should expect dispersion of these waves in the 
quantum gravity resulting from loop quantization techniques. 

It is well known that attempts of (loop) quantization of simplified models suffer from the problem 
that quantized reduced models have different features than the full theory - the more a system is 
simplified before quantization, the more likely that genuine features of full GR will be lost. In this 
paper we start, as a feasible compromise, from a slightly more general setting than pp-waves: The 



polarized Gowdy mode 
by Banerjee and Date 



in a form prepared for loop quantization, thoroughly studied and presented 
1CJ, Like space-times inhabited by pp-waves, this model is homogenous 
in two dimensions. But the Gowdy model has a different global topology. It is compact, whereas 
our model has the global topology of Minkowski space. Nonetheless due to homogeneity, we can 
choose an arbitrary finite area from the plane wave fronts and, as we consider only finite wave 
packets, the global topology of the inhomogeneous direction is not relevant for our purpose. 

The essential step in the reduction from the polarized Gowdy model to pp-waves is to single 
out space-times with waves traveling into one direction and so to avoid colliding plane waves 
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that lead to well-known complicated interaction processes and eventually to singularities |12[]. 
One such reduction of the polarized Gowdy model to pp- waves is carried out in [ijj. In that 
paper we used special coordinates introduced by Ehlers and Kundt [14]. This reduction involved 
second-class constraints and complicated Dirac bracket relations. Although the formulation is 
not a comfortable point of departure for loop quantization the Dirac brackets give analogs of the 

n 

commutation relation of a linear self-dual field [15(] with gravitational corrections that suggest 

corrections to h or the gravitational constant rather than a modified wave speed [13(] . In the 

present approach we formulate the reduction of the Gowdy model to pp-waves by means of a set of 

first-class constraints, so that the system becomes accessible to loop quantum techniques. A direct 

Hamiltonian formulation of pp-waves can be found in [171] . 

PI 

In this paper we review the formalism of Ref. [10] in the next section. Then we show in Section 
IIIII that the existence of a null Killing vector introduces constraints on the canonical variables. The 
new set of constraints for the model are shown to be first class in Section HVl Finally, initial steps 
toward quantization are carried out in Section [Vj 

II. THE POLARIZED GOWDY MODEL IN ASHTEKAR-BARBERO VARIABLES 

The vacuum Gowdy models are characterized by closed spatial topologies and two, spatial 

commuting Killing vectors. When the Killing vectors are orthogonal, then the model is "polarized". 

For our purposes it is convenient to employ the same formalism for pp-waves as used for polarized 

I I 

Gowdy models in Ashtekar-Barbero variables as formulated by Banerjee and Date [lOJ. The model 
has to be adapted to the setting of pp-waves. We replace the angular variable 9 on S 1 of the Gowdy 
model by the variable z on 3ft. We assume homogeneity in the x, y plane and wave propagation in 
the z direction. We choose adapted spatial triads with one leg in the z direction and two arbitrary 
orthogonal vectors in the homogeneous x, y plane. Densitized inverse triad variables are denoted 
by E a i, a = x,y and i = 1,2. As in [lo| we write 

£ := E z 3 (1) 
and introduce "polar coordinates" for the triad vectors in the x, y plane, 

E x i = E x cosn, E x 2 = E x smri, 

h (2) 

E y x = -Eysmr], EV 2 = EV cos 77. 

All variables depend only on z and the time variable t. The gauge group SU(2) of triad rotations 
has been reduced to U(l) rotations in the x,y plane, described by the dependence of the angle rj 



on z. In terms of these variables the spatial metric is 
ds 2 = h nh dx a dx b 



£ — dx 2 + £ — dy 2 + — — dz 2 . 
E x Ev £ 



(3) 
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Like the polarized Gowdy model of Ref . [101] the pp-wave space-times with either left of right 
moving waves are globally hyperbolic solutions to Einstein's equations that have two independent, 
orthogonal, commuting spatial Killing vectors. Unlike this Gowdy model, however, the pp-wave 
space-times are characterized by a null Killing field. To apply loop quantization techniques it is 
useful to introduce a hypersurface-orthogonal time coordinate and ask, under what conditions does 
the model describe pp-wave space-times? The answer will be in the form of a system of first class 
constraints in terms of the same phase space as the polarized Gowdy model. 

In the polarized Gowdy model the extrinsic curvature turns out to be diagonal, its compo- 
nents are denoted by K x , K y , and K z . The curvatures K x and K y turn out to be the corre- 
sponding diagonal elements of the Ashtekar-Barbero connection, divided by the Barbero-Immirzi 
parameter 7. As in Ref. [10], it is convenient to work with the 8-dimensional phase space 
{(K x , E x ), (K y ,E y ), (£,A), (77, P)}, where A is the component A z 3 of the connection, divided by 7, 
and the momentum P conjugate to r\ is constructed from the connection. The quantities E a , A and 
P transform as scalar densities, while the quantities K a , £ and 77 as scalars under diffeomorphisms 
along the z axis. The fundamental Poisson brackets are 

{K a (z),E b (z')} = K 5 b a 5(z-z'), {A(z),£(z')} = K 6(z-z'), {r,(z), P(z')} = KJ S(z - z') (4) 

where k = 87rGNewton is the gravitational constant. 

In terms of these variables we have the following set of first-class constraints for GR: The Gaufi 
constraint 
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G 



(5) 



(the prime denotes the derivative with respect to z), which generates rotations in the x,y plane, 
the diffeomorphism constraint 



K'E X + K'E y -£'A+-P 

x y 7 



(6) 



generating diffeomorphisms along the z-axis, and the Hamiltonian constraint 
1 



H 



K\/£E x Ey 



E x K x E y K y + (E X K X + EyR y )£ U+±)- \ £ ' 2 ~ ££ " 



E'y 

E x 



+ -££' '(In E x E y )' 



(7) 



W£E x Ey 



G 2 - 7 



£ 



E x Ey 



G 



In the following calculations it is convenient to use weakly equivalent forms of the last two 
constraints. Using the Gaufi constraint the diffeomorphism constraint can be rewritten in the form 



c = i 

K 



K'ji 1 ^K y E y -£' (A+^ 



+ V'G. 



(8) 



We make use of the weakly equivalent constraint 



C :=K' X E X + K' y E y -£' [A + 



kC. 



(9) 



This alternate form of the diffeomorphism constraint, C, is proportional to C modulo G. Similarly 
we define a form of the Hamiltonian constraint H, 
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H := —E x K x E y Ky — (E X K x +E y K y )£ ( A + £"+- 



lw2 
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weakly equal to k^J£E x EvH. 



III. REDUCTION TO PP- WAVES 



*«>.§); 



^££'{\nE x E y )' 
(10) 



We want to consider finite pulses of plane-fronted, parallel waves, traveling either in the positive 
or in the negative z-direction through flat space. Such waves are characterized by a null Killing 
vector field in the direction of propagation. In order to formulate the null Killing vector field we 
add an orthogonal time coordinate to the spatial manifold with the metric ([3]). Using a lapse 
function N = N(t, z) we have 



ds< 



E y 



E x 



E x E y 



-N 2 dt 2 + £ — dx 2 + £ — dy 2 + , 
E x Ev y £ 



■dz 1 



(11) 



The existence of the null Killing field satisfying V ' (^v) = gives rise to constraints on the phase 
space variables. Using the relations worked out in Appendix |A] we find two new constraints 



U x = E X K X --£'--£ 
2 2 

U v = E y K v --£' + -£ 
y y 2 2 



E yl 


E xl 


~E~y 


~~ "eF 


E yl 


E xl 


~Ev 


~ ~E* 



0, 
0. 



(12) 
(13) 



where the minus sign in k 11 in (|A2p was chosen [3]. The two relations render the diffeomorphism 
constraint and the Hamiltonian constraint equivalent; rewriting C and H in terms of U x and U y 
gives 

E y \'^ 



E x J 



^-(\uE x E y )' -£' (a + ^- 
2 V 7 



(14) 



The weak equality in (|14p implies that modulo U x , U v , and the Gaufi constraint 



E x Ev 



This means that if we choose the lapse function N = y e x ev ^en the time evolution generated 
by H is equivalent to a spatial diffeomorphism generated by C, i. e. time derivatives are equal to 
minus z-derivatives and the variables depend only on t — z. The waves travel without dispersion. 

In coordinates, such that g zz = 1, i. e. £ = E x E y and A = rj = 0, the constraint equation C = 
reduces to E x " E y + E x E y " = Q. In the formulation by Ehlers and Kundt 14] (see also Q), 



where E x = Le " and E y = Le^ , this becomes the the Einstein equation for pp-waves 

L" + L{f5') 2 = 0, (15) 
where the prime and the functional dependence are in terms of t ± z. 

IV. THE CONSTRAINT ALGEBRA 

We now investigate the possibility of imposing the pp-wave conditions (fT2j) and ([15]) as con- 
straints, augmenting the GR constraints G, C, and H. It turns out that the whole set of constraints 
is reducible and that the resulting reduced system is first class. 

The expressions U x and U y , smeared out with test functions, are 

U a [f] := J dzf(z)U a (z). (16) 
The Poisson bracket structure may be summarized in the following matrix 

WLfl.fsb]} {U,lf},U y {g}} J 2 \ -1 lJJ 
This matrix is diagonalized by the linear combinations 

U+ := U x + U y , U- := U x - U y , (18) 

explicitly 

U+ = E X K X + E y K y - £' and (19) 

/ E y \' 

V- = E X K X - E y K y Sim- . (20) 
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The algebra of these constraints is 



{[/+[/], U+[g}} = {U + [f],U4g}} = 0, and 
{[/_[/], U-\g]} = 2jdz(fg-fg')S. 



In local form, the non-vanishing Poisson brackets are 



{U„(z),U„(z r )} = -2[£(z)+£(z r )]6'(z-z'), 



(21) 
(22) 

(23) 



accordingly the constraints U-(z) are second-class. But the full set of constraints G, C, H, U + , and 
U- are reducible. In Appendix B we derive the relation 



U 2 _ = Ul + 28 



E X EV J V 7 



U+ + 4£U'. + AH- A£C 



(24) 



showing that U- is identically zero on the {C = 0, H = 0, = 0}-constraint surface. The 
constraint U- depends nonlinearly on these constraints. 

The constraint U- generates no further secondary constraints since the Poisson brackets of U- 
with G, C and H weakly or strongly vanish, 

1 



{U4f},G\g}} = 0, {U-[f},C[g}} = --U4f'g}*0, 

K 



{U4f],H[g]} = -U^ 

K 



£ f s+ £ ' /s 



E x Ev 



1 



U 4 



£ 



E x Ev 



fg 



(25) 



0; 



V£E x Ev 

U- is gauge invariant and conserved under spatial diffeomorphisms and time evolution 
Remarkably the constraint U+ weakly Poisson-commutes with G, C, and H, 

{U + [f},G[g]} = 0, 



{U + [f),C\g)} = — U + [f'g}^0, 

K 



{U + [f],H[g}} = -U^ 



£ 



E x Ev 



fg 



H[fg] « 0. 



(26) 
(27) 

(28) 



Thus, U+ can be added as another first-class constraint. The model now has the enlarged Poisson 
bracket algebra of the first-class constraints G,C, H, and £/+ that includes the standard algebra 
of constraints, adapted to the Gowdy model [10|, 



(29) 



{G[f],G[g}} = {G[f},H[g}} = 0, {G[f],C[g}} = -G[f'g], 



{C[f],C[g}} = C[fg' - fg], {C[f],H[g]} = H[fg'], 



{H[f],H[g]} = C 



(fg' - fg) 



E X EV 
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This sets the stage for quantization as now the algebra is first class, unlike in the analysis of Ref. 

0. 

The number of four first-class constraints is the maximum that can be imposed on a system with 
four canonical degrees of freedom. According to their nature as generators of gauge transformations 
the three sets of standard constraints of GR reduce the number of canonical degrees of freedom to 
one (two phase space functions). The constraints U+(z), on the other hand, reduce the two phase 
space degrees of freedom at every point - one field variable and its conjugate momentum - to one. 
By doing this £/+ reduces the physical degrees of freedom and but does not fall into the class of 
gauge generators, even though it is first-class, thus providing another counterexample to the Dirac 
conjecture. 

The quantity E x K x + E v K y appearing in U+, is the densitized trace of the extrinsic curvature of 
the wave fronts, when embedded into three-dimensional space-time. It has a geometric interpreta- 
tion as the expansion or contraction of an transverse area element under time evolution. Likewise, 
the difference E X K X — E y K y in the constraint f7_ has an interpretation as shear of transverse area 
elements. The analysis of the this section shows that area expansion can be integrated into the 
first-class constraint system, whereas shear cannot. 



V. PREPARATION FOR QUANTIZATION 

For quantization it is important that £/+ can be given a meaning as a well-defined operator. 
Indeed all the constraints can be formulated in a way that anticipates the construction of the 
corresponding loop quantum operators. Both E X K X + E y K y and £' are scalar densities, which 
can be naturally integrated along z, so to construct an operator we have to integrate them. The 
integral over some interval I of the coordinate z is 

U+ [1] = J dz(E x K x + E»K y ) -S+ + S-, (30) 

where £± are the values of £ at the endpoints of X. £ has a meaningful operator equivalent in the 

n 

adapted LQG framework of Ref. In analogy to full LQG the integral in (f30|) can be obtained 
as the Poisson bracket 

The first expression in the Poisson bracket is the first term of the kinetic term in the Hamiltonian 
constraint of (28) in Ref. The second expression is the volume of a sandwich of space, 
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constructed from a fiducial (unit) area in the x, y plane as basis and an interval I in the z- 
direction. Both expressions have operator equivalents in standard LQG, for the present case we 
find the corresponding operators in equations (31) and (32). Now we are in a position to 
express all first-class constraints in terms of loop quantum operators, acting on one-dimensional 



spin network states, as demonstrated in 



111- 



VI. CONCLUSION AND OUTLOOK 

In this paper we showed how the polarized Gowdy model can be further reduced to describe 
pp- waves. Somewhat surprisingly, the constraints from GR augmented by those derived from the 
Killing equations are reducible and that the reduced constraints form a first class system, which 
is amenable to loop quantum gravity quantization techniques. Much of the preparation of the 



quantization is carried out in Ref. U|. Work is underway on quantization of the model, which 
if all goes well, will provide a framework in which to perform ah initio investigations on Lorentz 
violation and deformation for pp-wave space-times. 

In classical GR the existence of a null Killing vector field in the direction of propagation of 
gravitational pp-waves guarantees dispersion-free propagation of such waves at a constant speed. 
In the present paper we show that a linear combination of two Killing equations is implemented 
as a first class constraint, which describes the expansion of null geodesies. The full information 
of the Killing equations, including shear, namely U+ and U-, that would render the Hamiltonian 
and the diffeomorphism constraint weakly equivalent, cannot be imposed in the form of first-class 
constraints. More specifically, if we exchange the first-class set G, C, H, U+ for G, C, £/+ and 
U-, to make H equivalent to C (see (fH|) ). the system becomes second-class. So, in order to be 
able to apply loop quantization techniques, we keep both C and H as independent constraints and 
the classical equivalence of time derivatives to space derivatives is not manifest in the constraint 
algebra. 

If we insist on implementing the full contents of the Killing equations as constraints, we are 
led to handle the second-class constraints U— with the method of Dirac brackets and impose them 
strongly. This program was carried out in detail in [jjj in a non-diffeomorphism invariant manner. 

An ongoing quantum investigation in the spirit of Ref. [11] will bring more clarity about possible 
dispersion of pp-waves in LQG. The results will be interesting to compare with Ref. [9J, which deals 
with gravitational wave dispersion in a cosmological background using perturbative techniques. In 
isolating pure gravitational wave propagation in one direction we have ended up with a more simple 
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model. With the higher degree of simplification we may expect complementary results to Ref. 



.a. 



Appendix A: Space-time metric, Killing equation 



In order to formulate a null Killing vector field we have to construct a space-time manifold. To 
this purpose we add an orthogonal time coordinate to the space-like manifold with the metric ([3]). 
The latter one becomes supplemented by a lapse function N = N(t, z), so that 



ds z = -N 2 dt 2 + £ ^- dx 2 + £ ^- dyV — - 
E x Ev £ 



E x 



E X E V 



■dz 1 



(Al) 



The non-vanishing components of the Levi-Civita connection are the following Christoffel symbols 



iV ' 



_ £NN> 

00 E x J^y > 



r o _ r o _N' 

1 03 — 1 30 — jy 7 ) 



33 



2N 2 V £ 



2E x Ev V £ 



2 V 5 



pO - 



1 f E y \ 

^t { £ E x ~) , r ° 1 = r i° = 2^ 



1 ^ * 



1 £ 



2E x Ev \ E x 



2 \ e x 



po - 

22 - 2Ar2 



1 £!/ 



1 £!/ / E x 
r 02 = r 20 = ~^=9t I £ 



2 5E a 



r 3 

1 22 



1 £ 



r 2 =l^L(e— V 

2E x Ev V Ev J ' ±23 32 2££ a; V / ' 



With the above four-metric a null vector field in the z-direction is of the form 

fc M = 



E x Ev 



£ 



fc,0,0,±iVfc 



(A2) 



with fc being a function of t and z. For fc^ to be a Killing vector field, the Killing equations 
^ uJ&v) = in terms of covariant derivatives with respect to the Levi-Civita connection must hold. 
From Vtfc* = we obtain the time derivative of fc, 

1 £ 



A; 

2E x Ey~"\ £ 
whereas from V z k z = we get the z-derivative, 

'E x Ey 



E x Ev N 



fc, 



(A3) 



1 



1 £ fE x Ey\' N' 
+ 2E x Ey\~£~) + lV 



2N V V s 

With these derivatives the equation V uk z -\ = reduces to an identity. 



fc. 



(A4) 
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Before imposing the remaining two non-trivial Killing equations V a /c a = we introduce the 
extrinsic curvature. In the case of a vanishing shift vector it is 

K a b = --^g'ab, (A5) 

so that we can express all time derivatives in the Christoffel symbols in terms of it. In the case 
of the polarized Gowdy model [l(J the diagonal components of the Ashtekar-Barbero connection, 
which are canonically conjugate to E x and E y , are - modulo the Barbero-Immirzi parameter - 
equal to the extrinsic curvature components. 

In the following we will use the components K a l = e^K^-, converted with the un-densitized 
radial triad components 



, P EV V 2 LEV rj 
e x = x j£—cos-, e x =^£— 8m - 



(A6) 



E* . rj 2 E* r, 3 E>E» 

t — — sm — , e v = \ £ — — cos - , e z = \ — - — , 
Ev 2 v V Ev 2 V £ 

corresponding to E x , E y , and £, respectively. In terms of these components we have 



1 E x „ / 'E"\ 



K x := y/(K x ^ + { K-Y = ^^d t [£- j (A7) 



Ky := ^/(V) 2 + (V) 2 = ^V^ ) m 

3 1 £ I E x E y \ 1 / E x Ev £\ 
2 " 2NE^ 1 \T~ ) ~2N\E^ + &~IJ- { > 

To express V a /c a = in terms of canonical variables we can use and T^ 2 m t ne form 



n K x £Ev n K v £E X /A . 

r?i = —— andrt 2 = S\— — . (A10) 
N V E x N V E y 



The Killing equations are independent of k 



T E X K X = l -£[ £ - + — - — ) (All) 
2 \£ Ev E x ' 



1 „ / '£' Ey' E 



which are the constraints (fT2l) and (fT3l), 
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Appendix B: The relation between £71, [/+, C and /J 



In this appendix we show some steps of the straightforward demonstration of the dependence 
of [/_, U + , C and -ff. From the expressions for f7_ and U + in (fT9]) and ([20]) we see that 

1 



and 



E X K V 



E y K,, 



E v \' 

U- + U-+S' + e [hx— J 



r - r - + ttV -^ ln — j 



(Bl) 



(B2) 



give 



(B3) 



and 



e*k z wk, = -u'i + - e'u + - j - 2 £ ) t 7 - 



m — 

E x 



Using these results in 77 produces 



H 



-8 8" + ^8 2 



1 



1 



£'E/+ + - E/"_ + - 8 In 



E x J 



4 2 

2 

~ 2 



#2/y 
IP 



L/_ - f M + — I 17+ 

7. 



i f £' (In E X E V )' -SS'(a+ 7 £\ 



Differentiating U+ and inserting into C gives 



1 1 / E y \' 

C = U' + -~QnE*E»yu + + ~(hi—) U. 



1 8' (In E X EV)' + 8" + ^8 



m tt" 



8'{A+^\. 
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(B4) 



(B5) 



(B6) 



From this we can eliminate 88' yA + 2- J from H. This substitution is sufficient to express U- 
completely in terms of the first-class constraints C, H, and E/+, resulting in relation (I24p . Thus, 
we have a quadratic equation for f/_, the solutions of which are of course also zero modulo the 
first-class constraints. 
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